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ABSTRACT 

We present results on the location, physical conditions, and geometry of the outflow in the Seyfert 
1 galaxy NGC 3783 from a study of the variable intrinsic UV absorption. Based on analysis of 18 
observations with the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope and 
6 observations with the Far Ultraviolet Spectroscopic Explorer obtained between 2000 February and 
2002 January, we obtain the following results: 1) The lowest-ionization species detected in each of 
the three strong kinematic components (components 1 - 3 at radial velocities —1350, —550, and 
—725 km s _1 , respectively) varied, with equivalent widths inversely correlated with the continuum 
flux. This indicates the ionization structure in the absorbers responded to changes in the photoion- 
izing flux, with variations occurring over the weekly timescales sampled by our observations. 2) A 
multi-component model of the line-of-sight absorption covering factors, which includes an unocculted 
narrow emission-line region (NLR) and separate covering factors derived for the broad line region and 
continuum emission sources, predicts saturation in several lines, consistent with the lack of observed 
variability in these lines. Differences in covering factors and kinematic structure imply component 1 is 
comprised of two physically distinct regions (la and lb). 3) We obtain column densities for the indi- 
vidual metastable levels from the resolved C III* A1175 absorption complex in component la. Based 
on our computed metastable level populations, the electron density of this absorber is ~ 3xl0 4 cm -3 . 
Combined with photoionization modeling results, this places component la at ~ 25 pc from the cen- 
tral source. 5) Using time-dependent calculations, we are able to reproduce the detailed variability 
observed in component 1 and derive upper limits on the distances for components 2 and 3 of < 25 
and < 50 pc, respectively. 6) The ionization parameters derived for the higher ionization UV ab- 
sorbers (components lb, 2, and 3 with log(t/)~ —0.5) are consistent with the modeling results for 
the lowest-ionization X-ray component, but with smaller total column density. The high-ionization 
UV components are found to have similar pressures as the three X-ray ionization components. These 
results are consistent with an inhomogeneous wind model for the outflow in NGC 3783, with denser, 
colder, lower-ionization regions embedded in more highly-ionized gas. 7) Based on the predicted 
emission-line luminosities, global covering factor constraints, and distances derived for the UV ab- 
sorbers, they may be identified with emission-line gas observed in the inner NLR of AGNs. We 
explore constraints for dynamical models of AGN outflows implied by these results. 
Subject headings: galaxies: individual (NGC 3783) — galaxies: active — galaxies: Seyfert — ultravi- 
olet: galaxies 
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1. INTRODUCTION 

Mass outflow, seen as blueshifted absorption in UV 
and X-ray spectra, is an important com ponent of active 
galac tic nuclei (see the recent review in lCrenshaw et aTl 
2003). This "intrinsic absorption" is ubiquitous in 
nearby AGNs, appearing in over half of Seyfert 1 galaxies 
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having high-quality UV s pectra obtained with t he Hubble 
Space Telescope (HST) l)Crenshaw e Lai] [1999) an d the 
Far Ultraviolet Spectroscopic Explorer (FUSE) QKriss I 
2002). Spectra from the Advanced Satellite for Cosmol- 
ogy and Astrophysics (ASCA ) identified X-ray "warm ab- 
sorbers" , modeled as absorption edges, in a similar per - 
centage of objects ijRevnolds Hl997tlGeorge et al.lll998ft . 
and many studies have explored the connection between 
the a bsorption observed i n the X-ray and UV bandpasses 
(e.g. iMathur et alJlT994j) . Large total ejected masses 
have been inferred for these outflows, exceeding the ac- 
cretion rate of the central black hole in some cases, in- 
dicating mass outflow may play a n important role in 
the overall ener getics in AGNs (e.g. IMathur et alJl!995l 
iRevnolds lll997|) . 

Variability in the intrinsic UV absorption in Seyfert 
galaxies is common. All objects with high-resolution UV 
spectra obtained at multiple epochs ex hibit substantial 
varia tions in their absorption strengths l|Crenshaw et all 
120031 and references therein). Absorption variability 
could result from: (a) a response to changes in the ioniz- 
ing AGN flux such that the total column density of the 
absorber remains constant but the ionization structure 
changes, (b) regions of condensation/evaporation in our 
line-of-sight to the AGN emission sources due to ther- 
mal perturbations, or (c) bulk motion of the absorber 
transverse to our line-of-sight. In the latter case, the 
observed equivalent widths could vary due to either a 
change in the line-of-sight covering factor of the back- 
ground emission, or a change in the total column den- 
sity of gas seen due to a shifting of different regions of 
the outflow across our sightline. In each case, the mea- 
sured variability characteristics provide important con- 
straints on the dynamics, geometry, and/or physical state 
of the AGN absorbers. For motion across the background 
AGN emission, the variability timescales can constrain 
the transverse component of the kinematics of the ab- 
sorber and the absorption-emission geometry. For ioniza- 
tion changes, the observed absorption variability can be 
used to constrain the gas number density and, combined 
with photoionization models, determine the distance of 
the absorber from the central source. These parameters 
are needed to determine the mechanism driving the mass 
outflow and the source of the absorption gas, and assess 
its overall role in the energetics of the AGN. 

The bright Seyfert 1 galaxy NGC 3783 has a rich UV 
and X-ray absorption spectrum. High-resolution obser- 
vations with the Goddard High Resolution Spectrograph 
(GHRS) and the Space Telescope Imaging Spectrograph 
(STIS) aboard HST showed the UV absorption is highly 
variable. Three distinct kinematic components of ab- 
sorption appeared independently over yearly timescales: 
components 1-3 having radial velocities v r w —1350, 
-550, and -725 km s" 1 and widths F WHM w 190, 170, 
and 280 km s _1 l|Kraemer et al.l2001l hereafter KC01). 16 
These long-term changes were found to be inconsistent 
with ionization changes, and thus interpreted as a sig- 
nature of transverse motion of the absorbers by KC01. 
High-resolution X-ray observations with the Chandra X- 
ray Observatory (CXO) revealed a spectrum with nu- 

16 Tentative detection of a w eak, fourth compon ent at v r £S 
— 1050 km s" 1 was described in Ga bel et alj 12003a!') : we do not 
treat this component in the following analysis. 



merou s absorption lines from a large range in ionization 
states l|Kaspi et al.ll200l|) . 

We have undertaken an intensive, multiwavclcngth 
campaign on NGC 3783 with ffST/STIS, FUSE, and 
CXO to monitor the absorption properties. Earlier pa- 
pers i n this series have presented studies of the mean 
X-ray <|Kaspi et al]l2002t Paper I) and UV ()Gabel et alJ 
l2003al Paper II) absorption spectra, analysis of a de- 
crease in radial vel ocity detected in UV component 1 
ijGabel et al.ll2003H Paper III), and varia bility and de- 
tailed modeling of the X-ray absorption l)Netzer et alJ 
I200.1l Paper IV). In this paper, we present a study of 
the variability and physical conditions in the UV absorp- 
tion based on analysis of the FUSE and STIS spectra. In 
§2, we review the observations and present the UV con- 
tinuum light curve; in §3, we present measurements of 
the absorption parameters and variability. We analyze 
metastable C III* A1175 absorption detected in compo- 
nent 1 in §4 to derive the number density in this absorber. 
In §5, we present detailed modeling of the UV absorp- 
tion, making use of observed variability in the spectrum, 
and explore time-dependent solutions in response to the 
continuum variations. We explore global models of the 
outflow in NGC 3783 in §6, based on the combined re- 
sults from the UV and X-ray analysis. 

2. OBSERVATIONS AND THE UV CONTINUUM LIGHT 
CURVE 

2.1. FUSE and HST/STIS Echelle Spectra 

We present results from a total of 18 medium- 
resolution STIS echelle spectra (SI - S18) and 6 FUSE 
spectra (Fl - F6) of the nucleus of NGC 3783, obtained 
between 2000 February and 2002 January. A detailed de- 
scription of the observations and data reduction is given 
in Paper II; here we present a brief overview. 

Each STIS observation was obtained using the 
0'.'2 x 0'.'2 aperture and the E140M grating, which spans 
1150-1730 A, and consisted of two HST orbits for a to- 
tal exposure time of ~ 4.9 ks (except SI which was 5.4 
ks). The STIS spectra were processed with IDL software 
developed at NASA's Goddard Space Flight Center for 
the Instrument Definition Team, which includes a proce- 
dure to remove background light f rom each order using 
a scattered light model devised bv iLindlerl |l999). Our 
measurements of the residual fluxes in the cores of satu- 
rated interstellar Galactic lines show the scattered light 
was accurately removed (see Paper II). The extracted 
STIS spectra are sampled in ~ 0.017 A bins, thereby 
preserving the full resolution of the STIS/E140M grat- 
ing (FWHM « 7 km s" 1 ). 

The FUSE spectra, covering 905 - 1187 A, were ob- 
tained through the 30" x 30" aperture. Each spectrum 
was processed with the standard calibration pipeline, 
CALFUSE. For each observation, the eight individual 
spectra obtained with FUSE from the combination of 
four mirror/grating channels and two detectors, were 
coadded for all exposures. We corrected for nonlinear 
shifts in wavelength scale between the spectra from dif- 
ferent detector segments by cross-correlating over small 
bandpasses. The absolute wavelength scale was deter- 
mined by matching the velocities of Galactic lines in 
the FUSE spectrum with those in the STIS bandpass. 
Mean residual fluxes measured in the cores of saturated 
Galactic lines are consistent with zero within the noise 
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(i.e., standard deviation of the fluxes) in the troughs of 
these lines, indicating accurate background removal for 
all observations except F5. For this observation, we fit 
the remaining residual background flux in strong inter- 
stellar lines to match the other epochs, and subtracted 
the fit. The spectra were resampled into ~ 0.02 - 0.03 
A bins to increase the signal- to- noise ratio (S/N) while 
preserving the full resolution of FUSE, which is nomi- 
nally FWHM w 20 km s _1 . 

2.2. The UV Continuum Light Curve 

To put the continuum variations of our observations 
in perspective, Figure la shows the continuum flux light 
curve at 1470 A for all UV spectra of NGC 3783 obtained 
over the last 22 years. These observations were obtained 
with the International Ultraviolet Explorer (IUE) and 
the Faint Object Spectrograph (FOS), GHRS, and STIS 
(in low dispersion) on HST, which are listed in Table 1. 
We obtained the most recently processed version of each 
spectrum from the Multimission Archive at the Space 
Telescope Science Institute. We measured the contin- 
uum fluxes in each spectrum by averaging the elements 
in a 30 A bin centered at 1470 A in the observed frame, 
which is free of contamination from line emission. The 
lcr flux uncertainties were determined from the standard 
deviations. For the IUE spec tra, this technique is known 
to overestimate the errors (|Clavel et all I1991D: there- 
fore, f ollowing the procedure described in Kraemer et alJ 
( 2002), we scaled these uncertainties by a factor of 0.5 to 
ensure that observations taken on the same day agreed 
to within the errors on average. Due to the small wave- 
length coverage of the GHRS spectra, no continuum re- 
gions free of broad line emission were observed. To esti- 
mate the GHRS continuum fluxes, we used separate fits 
to the continua and broad emission lines from the STIS 
spectra and tested different linear combinations of these 
fits until we obtained an accurate match to the observed 
GHRS profiles. Figure la shows the STIS echelle obser- 
vations (JD w 2,452,000) sampled the UV continuum in 
a range of moderately-low to moderately-high flux states 
compared to the long-term light curve. 

In Figure lb, the light curve for the monitoring cam- 
paign observations is shown in more detail. We have in- 
cluded estimated fluxes from the four FUSE observations 
that were not simultaneous with STIS observations (Fl, 
F3, F5, F6) by extrapolating from the far-UV bandpass 
to 1470 A. The first four observations (SI - S4) were ob- 
tained at intervals of several months, the intensive mon- 
itoring phase (S5 - SI 7) sampled the continuum at 3 - 
8 day intervals, and the final observation (S18) was ob- 
tained about 9 months later. The extrema in continuum 
flux in the STIS observations occurred during the long- 
term sampling; S2 and S3 observed the continuum in the 
highest flux state and the final observation found it in the 
lowest state, with a peak amplitude variation of a factor 
of ^2.5. We identify two general phases during the in- 
tensive monitoring: the first four observations (S5 - S8) 
sampled the continuum in a relatively low-state, after 
which it increased by up to a factor of 1.7 (S10 and S12) 
over a period of about 2 weeks. Three FUSE observa- 
tions (F2 - F4) were obtained during the low-state STIS 
observations, while observation F5 observed the contin- 
uum in the highest state (four days after S12). In the 
subsequent variability analysis, we compare observations 



between mean low and high states derived by averaging 
representative spectra in each state. 

3. MEASUREMENTS OF ABSORPTION PARAMETERS 

A key issue in the analysis of AGN outflows is extract- 
ing accurate ionic column densities from the observed 
absorption lines, since they provide the basis for de- 
termining the physical conditions in the absorbers (i.e., 
the ionization state, total gas column, number density). 
The UV absorbers typically only partially occult the 
background AGN emission, as determined by the rel- 
ative strengt hs of the individual members of absorp- 
tion doublets llWamnler et al.lll993HHamann et al.lll997t 
iBarlow fc Sargent Ill997j) . and the effects of covering fac- 
tor (C) and optical depth (r) must be deconvolved to 
determine the column densities. In some cases, the C — t 
solution for observed features can be complicated due to 
blending of physically distinct absorption components or 
complex coverage of the background emission. Thus, it 
is useful to first explore absorption variability by com- 
paring equivalent widths. 

3.1. Variability in Absorption Equivalent Widths 

Figures 2a - c show the equivalent widths for key lines 
in the STIS spectra of components 1-3, respectively, 
plotted as a function of the UV continuum flux. The 
error bars represent our estimated measurement uncer- 
tainties, which are a combination of uncertainties due to 
spectral noise and fitting the intrinsic (i.e., unabsorbed) 
fluxes, via propagation of errors. We estimated uncer- 
tainties in the intrinsic fluxes by testing different empir- 
ical fits over the absorption features and selecting the 
range of what we deemed to be reasonable line profile 
shapes. We note that these derived uncertainties may in 
some cases be too conservative, as indicated by compar- 
ing the scatter in measured values with the error bars 
plotted in Figure 2 (e.g., N V in components 1 and 3). 
This may be due to systematic trends in how the intrin- 
sic flux was fit in each observation and/or overestimates 
of the uncertainties associated with the fitting. 

In each kinematic component, the lowest-ionization 
species with detectable absorption other than H I (Si IV 
in component 1; C IV and N V in component 2; 
C IV in component 3) show clear variations that are 
inversely proportional to the continuum flux. This is 
just as expected for unsaturated lines from relatively low- 
ionization species in a photoionized g their column 
densities decrease in higher flux states due to the in- 
creased ionization. This gives direct evidence that the 
ionization structure in these absorbers is dominated by 
photoionization from the central source. It also indi- 
cates these lines are not highly saturated, at least in the 
epochs with weaker observed absorption. To test the 
dependence of absorption variability on continuum flux 
more rigorously, we did a linear regression fit for each 
line. Figure 2 shows the results of the fits and gives the 
ratio of the computed slope to the la uncertainty in the 
slope, m/Am, to characterize the correlation. Fits to all 
of the lines listed above have a non-zero slope at greater 
than a 3cr level. Since the absorption strengths will not 
necessarily vary linearly with continuum flux, we also 
tested the Spearman rank correlation. All of the above 
lines show correlation with the flux at high significance, 
with probabilities of no correlation computed to be be- 
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tween 0.01 - 0.3 % from this test (values are listed in 
Figure 2). 

Conversely, the absorption strengths of C IV in com- 
ponent 1 and N V in components 1 and 3 did not vary 
strongly during the monitoring. This could indicate ei- 
ther these lines were saturated, with partial covering fac- 
tor, or their ionic column densities were not varying. In 
the latter case, the lack of variability could be due to ei- 
ther a relatively low electron density such that the ionic 
populations are unable to respond to continuum changes, 
or to the ionic structure of the gas, with these ions near 
the peak ionization states of their parent elements. These 
possibilities are explored below. We note that since each 
component has some lines that exhibit no variability (see 
Figure 2 and §3.2), the lines that do vary must result 
from changes in ionic column densities rather than cov- 
ering factors. 

3.2. Effective Line-of-Sight Covering Factors and 
Column Densities 

The doublet method for measuring the absorption pa- 
rameters provides a solution to a single C and r (at 
each resolution element in an absorption profile). How- 
ever, since the background AGN emission is comprised 
of multiple, physically distinct components with differ- 
ent sizes and geometries (i.e., a featureless continuum 
source and multiple kinematic components of line emis- 
sion), additional constraints may be needed to deter- 
mine the required effective covering factors, which are 
weighted combinations of covering factors of the differ- 
ent emission sources (Ganguly et al. 1999). Implicit in 
the doublet solution is that all emission sources have the 
same covering factor. In Paper II, we used the Lyman 
series lines to separate the continuum and emission line 
covering factors for the H I absorption. Here, we use 
variability in the background AGN emission as an addi- 
tional constraint to explore the effect of the NLR on the 
derived covering factors and obtain a consistent model 
to measure absorption column densities for all lines. 

3.2.1. Isolating the Narrow-Line Region Emission 

Profile 

We first compare spectra in different flux states to iso- 
late distinct emission-line components based on variabil- 
ity in the overall profiles. Figure 3a shows the C IV 
line profiles in mean high-state (S2 and S3) and low- 
state (S5 - S8) spectra; these observations were se- 
lected for comparison because they show the largest dif- 
ference in emission-line flux. The continuum flux has 
been subtracted from each spectrum and the low-state 
profile has been scaled by a factor of 1.4 to match the 
flux in the high-velocity wings of the high-state spec- 
trum. The profiles match very well at radial velocities 
\v r \ >1500 km s _1 but diverge at lower velocities, with 
increasing discrepancy toward line center. This effect is 
consistent with the superposition of a varying broad com- 
ponent and a non- varying narrower emission-line compo- 
nent, hereafter the BLR and NLR respectively. 

Assuming no change in the NLR flux between obser- 
vations, we isolate its profile by solving the following ex- 
pression for Fjvl_r at each radial velocity: 

Fi ow (v) x f sc - F high (v) = F NLR (v) x (f sc - 1), (1) 
where Fi ow and F^igh are the total observed emission- 
line fluxes in each state (i.e., the NLR + BLR fluxes) 



and f sc is the scale factor equating the low-state BLR 
flux with the high-state. The resulting NLR profile is 
plotted as a dotted line in Figure 3a. This analysis as- 
sumes the BLR scales by a uniform factor at all radial 
velocities between states. Equation 1 is not valid in spec- 
tral regions affected by absorption. We have fit the NLR 
profile using Gaussians for each of the C IV doublet lines 
constrained to have their intrinsic 2:1 flux ratio and ve- 
locity separation, giving a width of a = 500 ±120 km s _1 
{FWHM = 1180 ±280 km s" 1 ) at a radial velocity of 
±50 ±60 km s _1 with respect to the systemic velocity. 
The combined fit, which appears symmetrical due to the 
broad line width relative to the separation of the doublet 
members, is plotted over the residual NLR flux in Figure 
3a. 

We applied the same analysis to derive the NLR pro- 
files of other emission lines using equation 1. Results 
for Lya, N V, and Si IV (f sc =1.4, 1.8, and 1.8, respec- 
tively) are shown in Figure 3b. A scaled template of 
the Gaussian fit to the C IV NLR profile is overlaid on 
each residual NLR profile, including both lines for the 
N V and Si IV doublets. The discrepancy in the residual 
fluxes in the blue wing of N V (v r < —2000 km s _1 ) and 
the red wing of Lya (v r > 3000 km s _1 ) is due to the 
different BLR flux scale factors for the two lines and the 
overlap in their BLR emission. The excess emission in 
the Si IV profile between ~ 1000 - 2500 km s" 1 is due 
to the O IV] emission-line multiplet. The Lya NLR is 
not well matched by the C IV template, appearing nar- 
rower in the unabsorbed red wing of the profile. Figure 3 
shows the derived NLR fluxes contribute substantially at 
the wavelengths of many of the absorption features and 
thus its covering factor could affect the measurement of 
column densities, which we explore below. 

There are some caveats regarding the derivation of the 
NLR profiles. First, the assumption that the BLR scale 
factor between the low and high states is independent of 
radial velocity may introduce an error into the solution. 
Intensive UV - optical monitoring of continuum and BLR 
variability in the Seyfert 1 galaxy NGC 5548 has re- 
vealed its BL R variability is consistent with a constant 
virial product ({Peterson fc Wandel ll999UPeterson et alJ 
l2004|) . A consequence of this is that the integrated 
BLR profile becomes narrower in higher flux states; if 
this is the case for the UV lines in NGC 3783, then a 
velocity-dependent scale factor in equation 1 would be 
required to fully separate the NLR and BLR. Addition- 
ally, the profiles derived from this analysis are broader 
than are typical for UV NLR lines in Seyfert 1 galax- 
ies. For example, HST observations of NGC 5548 and 
NGC 4151 while in low flux states with little contami- 
nation from BLR emission revealed FWHM ss 500 and 
300 km s -1 , respectively, for the C IV NLR. However, 
the width derived for NGC 3783 is comparable to the 
rather broad NLR features in the Seyfert 2 galaxy NGC 
1068 llDietrich fc Wagner lll998t iKraemer et al.lll998aj) . 
To assess the NLR fit, we compare their flux ratios rel- 
ative to the measured [O III] A5007 line with those from 
other AGNs. From optical spectr a of NGC 378 3 ob- 
tained with a 5"xl5" aperture bv lEvansI l)1988f) . the 
ratio of C IV from our fit to the total [O III] flux is 
~0.5. In comparison, measurements of NGC 5548 yield 
a C IV : [O III] ratio of w 1.2, o btained with I" (C IV ) 
and 4"xl0" ([O III]) apertures l|Kraemer et al] ll998b). 
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Based on this, the C IV flux from our derived profile 
is reasonable. In contr ast, in a sample o f more lumi- 
nous radio-loud AGNs, iWills et all l)1993|) detected no 
UV NLR lines, using the observed [O III] line as a tem- 
plate, with upper limits on the C IV : [O III] flux ratios 
of w 0.1 -0.5. 

3.2.2. Covering Factor Model 

For absorption features imprinted on multiple discrete 
background emission sources, the normalized flux for the 
j th line can be written: 

I, = URtttfe-T' + 1 - C% (2) 

where the \ th individual emission source has fractional 
contribution to the total intrinsic (i.e., unabsorbed) flux, 
R) = F-/T, t [F-], and line-of-sight covering factor, Cj 
l|Gabel et al.l l2005'). The effective covering factor associ- 
ated with each line is the weighted combination of the 
individual covering factors: 

C^ZiiRiq}. (3) 

T hese equations a r e exte nsions of the expressions given 
in iGangulv et al. (1999) for the continuum and BLR 
to include an arbitrary number of background emis- 
sion sources. Combining equations 2 and 3 and solv- 
ing for r gives the fa miliar expression for optical depth 
l|Hamann et al.lll997j) . 

^T^TTc-)- (4) 

The total column densities for each line are then obtained 
by integrating 

N ^ v) -^Tx T ^ (5) 

over radial velocity ijSavage fc Sembachl I1991|) . where 
the optical depths are derived in each velocity bin from 
equations 3 and 4. 

Analysis of the Lyman lines in Paper II revealed the 
emission-lines are only partially occulted by the UV ab- 
sorbers. Since the NLR is generally much more extended 
than the BLR and continuum source in AGNs, we as- 
sume here a 3-component geometrical model in which 
the NLR is entirely unocculted by all components of UV 
absorption. Later, we explore this assumption based on 
geometrical constraints from our analysis. In Paper II, 
the emission-line and continuum covering factors were 
separated assuming a single emission-line region. Here, 
we incorporate the NLR - BLR separation for our 3- 
component coverage model (hereafter 3-C) by solving for 
the continuum and BLR covering factor profiles, C c and 
C BLR , using the Lyman lines as in Paper II after first 
subtracting a model of the Lya NLR. Given the mis- 
match between the Lya NLR residual and C IV template 
(see discussion above and Figure 3b) , we fit the Lya NLR 
with a narrower Gaussian that matched both the uncon- 
taminated red wing and core, and was constrained to be 
below the residual flux in the component 2 absorption 
feature (a =350 km s _1 ). Given the uncertainty in the 
NLR profile, we cannot rule out that the deep compo- 
nent 2 feature partially occults the NLR of Lya. 

Using equation 3, we calculated effective covering fac- 
tors for each absorption line, using the derived C c and 



c blr pron i eS; anc i (jNLR = q at all racua i velocities. Fig- 
ure 4 shows the resulting normalized unocculted flux lev- 
els (1 — C) for this 3-C model compared to the observed 
normalized absorption profiles. Regions where the resid- 
ual fluxes in the absorption troughs are close to the un- 
occulted flux levels (I 1 — C) indicate the lines are near 
saturation for this model. For each line, results for both 
low and high-state spectra are shown to demonstrate the 
variability and implications of the covering factor on ob- 
served line strengths. To increase the S/N in each spec- 
trum, we averaged multiple observations for each state, 
selecting observations that exhibited the largest varia- 
tions in absorption strengths and continuum flux during 
the intensive monitoring. For the STIS lines, the low 
state shown is the mean of S5 - S8, and the high state 
S10 and S12, except for Si IV, which had the weakest ab- 
sorption in observations S15, S16. For lines in the FUSE 
spectrum, the low-state is the mean of F2 - F4, and the 
high-state F5 - F6. Column densities for the low and 
high-states in each component computed with equations 
4 and 5 are given in Table 2. Below, we describe the 
absorption lines for each kinematic component. 

For the N V absorption, both doublet members are free 
of contamination over most of the absorption profiles, 
providing a comparison between the doublet solution to 
the covering factor (Cd) and our 3-C model to test the 
effect of an unocculted NLR. The unocculted flux levels 
for the doublet solution (1— Cd) are shown in Figure 4a 
with horizontal tick marks for N V. These values were 
derived in the cores of each kinematic component in the 
merged spectrum (Paper II). Close comparison reveals 
important differences in the implied column densities be- 
tween the two covering factor models. In components 1 
and 3, the doublet solution gives unsaturated absorption. 
This is seen in Figure 4a, where I > 1 — Cd in the weaker 
doublet member (N V A1242). The resulting N V col- 
umn densities are ~8 x 10 14 cm~ 2 (component 1) and 
1.5 x 10 15 cm -2 (component 3), with equal values in 
the low and high-state spectra within measurement un- 
certainties. Thus, the doublet solution implies the N V 
absorption in components 1 and 3 is unsaturated, with 
the above column densities, and without any detectable 
variability between flux states. A similar result is found 
consistently for all STIS observations. In contrast, the 
3-C coverage model implies N V is near saturation in 
these components, in both flux states, since Z w 1 — C 
consistently for both doublet lines. The primary differ- 
ence is due to the NLR flux underlying the red doublet 
members of these components (see Figure 3b). The lack 
of substantial variability in the equivalent widths (Figure 
2a and 2c) provides independent evidence for saturation 
of N V in these components. N V in component 2 is 
unsaturated for both covering factor models, consistent 
with the variability detected in its equivalent width. 

3.2.3. Component 1 

Component 1, at v r = —1350 km s , exhibits a rich 
absorption spectrum. In addition to the common C IV, 
N V, and O VI doublets, it includes lines from the rela- 
tively low- ionization species Si IV and C II, P V All 18 
(the P V A1128 doublet line is contaminated with Galac- 
tic absorption and unmeasurable) and the metastable 
C III* A1175 complex (Paper II). None of these lines 
is detectable in the other kinematic components in NGC 
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3783, and they are only seen in a small fraction of in- 
trinsic absorbers in AGNs (NGC 4151 is the only other 
Seyfert with reported C III* and P V absorption, Bro- 
mage et al. 1985, Espey et al. 1998; Si IV appears in ~ 
40% of Seyfert absorbers, Crenshaw et al. 1999). The 
C III* A1175 absorption gives a tight constraint on the 
number density of the absorber, as shown in §4. H I is 
detected in the Lyman series up to Ly#, while contam- 
ination from Galactic absorption prevents detection of 
the higher order lines; thus the H I column density listed 
in Table 2 should be considered a lower limit. 17 

There is evidence that multiple, physically distinct ab- 
sorption regions are overlapping in kinematic compo- 
nent 1 (KC01; Paper II). In Paper II, it was shown the 
Si IV covering factor from the doublet solution is lower 
than that derived from the Lyman lines and the N V dou- 
blet. This is confirmed for the new 3-C model: the col- 
umn density measured for the red Si IV doublet line with 
this covering factor is three times greater than the blue 
line, indicating the actual Si IV covering factor is signifi- 
cantly lower. In contrast, the 3-C model gives consistent 
results for the two N V doublet members: saturation in 
the blue- wing and core in both flux states (see above dis- 
cussion), and similar column densities in the red-wing 
of the profile where the absorption depths diverge from 
the unocculted flux levels (i.e., the absorption is unsatu- 
rated). Similarly, C IV A1548 and O VI A1038 are consis- 
tent with being saturated in their blue wings and cores 
with the 3-C model (the other doublet members for these 
ions are unmeasurable due to contamination with other 
absorption), and C IV diverges substantially in the red 
wing (O VI A1038 is contaminated with a detector arti- 
fact at these velocities and cannot be tested). Addition- 
ally, there are ion-dependent structural differences in the 
absorption profiles. As described in Paper II, the red 
wings of some lines, particularly N V and Lya, extend 
to lower velocities than Si IV and P V. Interestingly, it is 
at these velocities that the N V and C IV profiles diverge 
from the unocculted flux levels from the covering factor 
model. 

Thus, we treat the absorption from this kinematic re- 
gion as coming from two physical components: com- 
ponent la has relatively low covering factor and low- 
ionization and gives Si IV, C II, C III*, and P V; compo- 
nent lb is more highly ionized and has higher covering 
factor, contributing strongly to C IV, N V, O VI and 
the Lyman lines. The latter lines will have contribu- 
tions from both components, with absorption from com- 
ponent l a buried in the higher covering factor absorber 
(e.g., see iKraemer et all l2003|) . To measure the Si IV 
column density, we used the covering factor derived from 
the doublet pair (C =0.35; Paper II). For the other lines 
associated with component la, there is no independent 
measure of the covering factors (it is not possible to sep- 
arate the emission-line and continuum covering factors 
for this component from the single Si IV doublet). We 
assumed C =1 for these lines because their underlying 
emission is predominantly continuum flux. Since they 
are all weak (see Figure 4), the measurements will not 
be too far off unless C c is very small in this component. 

17 In Paper II, we missed detection of hj8 (as well as Ly£) in 
component 1; thus the value quoted there, which was measured 
from the Lye line, is smaller than the present study. 



For the O VI, N V, and C IV lines, we adopted the 3- 
C model. Since all these lines show saturation in this 
component over much of the profiles, we derived lower 
limits on their integrated column densities by adding es- 
timated uncertainties to the normalized fluxes, / + Al, 
in equation 4. Measured column densities and limits for 
the low and high-states are listed in Table 2. For mod- 
eling purposes (§5.2), we also measured the C IV and 
N V absorption in the red wing of the profile, where the 
lines are unsaturated and there is no contribution from 
component la based on the Si IV profile (v r — —1270 - 
-1170 km s" 1 ). This gives 3. 7xl0 14 and 1.0xl0 14 cm~ 2 
for the N V and C IV column densities, respectively. 

As described in §3.1, the Si IV absorption varied, with 
a general inverse correlation with continuum flux. With 
the relatively low Si IV covering factor, the moderate 
equivalent width variations translate into peak changes of 
a factor of ~ 4 between the mean low and high states dur- 
ing the intensive phase defined above. The column den- 
sity decreased somewhat less in the high-state S10 and 
S12 observations relative to the mean low-state (a factor 
~ 2). Figure 4c shows the same trend for C II: weak 
absorption is detectable in the ground state (C II A1335) 
and fine-structure (C II* A1336) lines in the low-state 
spectrum, but was not distinguishable from noise in the 
high-state. The weak P V A1118 and C III* A1175 ab- 
sorption appears to decrease as well, with only ~1.5cr 
detections in the high-state, though the limited S/N in 
the spectra of these lines makes this less certain. No sig- 
nificant variability was detected in Ly#, nor any other 
Lyman series lines, between flux states. 

3.2.4. Component 2 

In component 2 (v r = —550 km s _1 ), the N V and 
C IV absorption is relatively weak and variable (Figures 
2b and 4a). Our adopted covering factor model gives 
unsaturated absorption in these doublets, consistent with 
the observed variability. In contrast, Figure 4b shows 
the O VI A1038 line is strong and matches the model 
unocculted flux levels nearly identically across the entire 
profile in both states, consistent with heavy saturation of 
this line. Lyman line absorption is measurable only up to 
Ly7 in component 2 due to strong contamination of the 
higher order lines with Galactic absorption. The Ly7 line 
does not vary significantly, thus we adopt the H I column 
density measured from this line as a lower limit. Figure 
4b shows both the C III A977 and N III* A991 lines are 
contaminated with moderate Galactic H2 absorption at 
component 2 velocities, with neither line showing strong 
evidence for the presence of intrinsic absorption. Upper 
limits on these lines measured after the removal of the 
H2 model are given in Table 2. 

3.2.5. Component 3 

In component 3 (v r = —725 km s _1 ), the primary con- 
straints for modeling are the C III A977 line and the 
C IV doublet. Both lines exhibit variability, decreasing in 
epochs with higher continuum flux. The C III line is rela- 
tively strong in the low-state FUSE spectrum, but is sim- 
ilar to the noise level in the high state (Figure 4b), giving 
only an upper limit. Figure 4a shows that for the adopted 
covering factor model, the C IV A1551 line is near sat- 
uration in the low-state, but well above the unocculted 
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flux level in the high-state. Thus, for this model, the 
equivalent widths in C IV vary only moderately (Figure 
2c), but the variation in column density is large between 
flux states due to the low effective covering factor (C « 
0.35 in the core of component 3). The C IV AI548 line 
is blended with the component I C IV A1551 line, thus 
no comparison of our adopted model with the doublet 
solution is possible. The O VI A1038 and N V doublet 
lines (see above) are consistent with being saturated in 
both flux states, based on the unocculted flux levels. The 
highest order Lyman line not contaminated with other 
absorption in component 3 is Lye. Since it does not vary 
significantly between flux states, we adopt the resulting 
H I column density from this line as a lower limit. The 
N III A989 component 3 line coincides with the strong 
N III* A991 component 1 line, and cannot be measured. 
N III* A991 in component 3 is moderately contaminated 
with Galactic H2. Given the limited S/N in this region, 
we take the measurement of N III* after the removal of 
the H2 model as an upper limit. 

4. CONSTRAINTS ON THE DENSITY FROM THE 
METASTABLE C III ABSORPTION 

The C III* A1175 multiplet lines have been used as a 
density diagnostic for AGN absorbers in sever al studies. 
However, as pointed out bv lBehar et al.1 l)2003|) . the high 
densities derived in these studies, including our Paper II, 
were based on calculations of level populations that only 
treated the 3 Pi level. The J =0 and 2 levels have much 
lower radiative transition probabilities to the ground 
state (A «0 and 5xl0 -3 s _1 , respectively) than the J =1 
level (A =75 s _1 ), and thus are populated at densities 
that are lower by several orde r s of magnitude (see Figur e 
1 in lBhatia fc Kastner lli993HKastner fc Bhatia l [l9921. 
This has important consequences for the interpretation 
of the absorbers since the gas density is needed to de- 
termine the location and physical depth of the gas. To 
address this, we computed the relative populations of the 
3 Pj levels for a range of densities and electron temper- 
atures, T e , that are expected for the photoionizcd AGN 
outflows seen in the UV. Our calculations extend the 
results of Bhatia & Kastner (1993), who give results ap- 
plicable to collisionally ionized plasmas with T e > 40000 
K, to lower temperatures. Collisional excitation and de- 
excitation and radiative decay between all levels were 
treated. We treated the 6 lowest terms/levels of the C +2 
ion: the ground state 2s 2 1 So term, the three 2s2p 3 P 
levels, 2s2p 1 Pi, and 2p 2 3 P. Temperature dependen t col- 
lision strengths from lBerrington et al.1 1)1985111989^ were 
used for transitions between levels and terms, respec- 
tively, where we have interpolated between their listed 
temperatures. Radiative transition rates were obtained 
from lBhatia fc Kastner I l)1993j) . lMorton I l)1991[K and the 
NIST Atomic Spectra Database website. 

Figure 5 (top panel) shows the computed populations 
over a large range in density for T e = 16000, 20000, and 
40000 K. This shows that the C III* A1175 absorption 
complex can serve as a powerful probe of the physical 
conditions in the absorber. The relative populations of 
the three levels are very sensitive to the electron density, 
but insensitive to temperature. This is seen in the bot- 
tom panel of Figure 5, where the ratio of the J =2 : J =0 
level populations are plotted for the temperatures shown 
in the top panel. Additionally, the absolute populations 



of the 3 P levels are very sensitive to the gas temperature. 

In practice, full utilization of these lines as diagnos- 
tics requires that the absorption features are sufficiently 
narrow so that the individual lines in the complex are 
not heavily blended, and sufficient resolving power to 
separate the lines. These conditions are met with the 
STIS spectra of component 1 in NGC 3783 as seen in 
Figure 6, which shows only mild blending of the indi- 
vidual lines. The location of the six multiplet lines are 
marked and identified by the J level of the transition. 
A fit to the C III* A1175 complex is shown as a dashed 
line, using the width and centroid of the Si IV absorp- 
tion profile; best-fit column densities for each level are 
given below the spectrum. The column density ratio of 
the J =2 : J =0 levels, Nj= 2 /Nj= = 2.9±1.4, gives 
n e =3~ti 5 X 10 4 cm~ 3 , independent of temperature. The 
lack of detection of the A1174.93 line is consistent with 
this density since Figure 5 shows the population of the 
J = 1 level will be approximately four orders of magni- 
tude lower than the other levels. Use of these level popu- 
lations as a temperature diagnostic requires knowing the 
total abundance of the C +2 ion; since the C III A977 line 
is unmeasurable due to contamination (and saturated for 
the implied column density), this requires the results of 
photoionization modeling, which is presented in §5. 

Our calculations do not treat the effects of the radia- 
tion field on the level/term populations. This could de- 
crease the ground state or metastable level populations 
due to continuum pumping, or increase the relative pop- 
ulations of the metastable levels due to recombination 
followed by cascade to the 2s2p levels. We have com- 
puted photoionization models with Cloudy to test this 
and find it has a negligible effect on the metastable level 
populations in this component. 

5. PHOTOIONIZATION MODELS 
5.1. Input Parameters and Assumptions 

We compare the observed ionic column densities with 
predictio ns from the photoionization modeling code 
Cloudy l|Ferland et al.lll998t) to constrain the physical 
conditions and location of the intrinsic absorbers. The 
absorbers are assumed to be uniform plane parallel slabs 
of constant density, photoionized by the AGN at a dis- 
tance R from the central source. These calculations ap- 
ply for a gas in ionization and thermal equilibrium with 
the ionizing flux; this is explored in §5.3 based on the 
observed variability in the absorption. The models are 
specified by the spectral energy distribution (SED) of the 
ionizing continuum, the total hydrogen column density 
(Nh) and elemental abundances in the absorber, and the 
ionization parameter (U = Q /AirR^nHc), which gives the 
ratio of the density of H ionizing photons at the face of 
the absorber to the gas density, nu- 

The choice of SED for the unobservable ionizing con- 
tinuum is a source of modeling uncerta inty in AGNs (e.g. 
iMathur et aT1ll994 iKaspi et alJl2001^ . To facilitate the 
comparison with the earlier UV results, we adopted an 
SED based on the one used by KC01. This consists of 
multiple power-law components {F v cx v~ a ) constrained 
by the UV and X-ray observations. The UV power-law 
component is based on continuum-flux measurements in 
the two epochs with simultaneous STIS and FUSE ob- 
servations. After correcting the spectrum for Galactic 
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extinctio n (E(B - V) = 0.1 19; iReichert et all fl99^) 
with the iCardelli et al.1 1(19891) reddening curve, the con- 
tinuum flux ratio measured at 950 A and 1690 A gives 
a spectral index of a^ v ~ 1. As discussed in KC01, ex- 
trapolating the UV power-law to X-ray energies far over- 
estimates the observed flux, thereby requiring a spectral 
break in the unobservable EUV - soft X-ray flux to con- 
nect the absorption in the two bands. Thus, we extrap- 
olated from the observed UV flux at the Lyman limit to 
the observed flux at 0.6 keV in the Chandra spectrum, 
giving a spectral index af uv — 1.4. For the hard X-ray 
index, we used the value derived from the merged Chan- 
dra spectrum in Paper I, a„ rs 0.7 (T =1.7). For this 
SED and the adopted d istance of 39 Mpc for NGC 3783 
(using z =0.00976 from lde Vaucouleurs et al.lll99lL and 
H a =75 km s _1 ), this gives an H-ionizing photon lu- 
minosity Q =1.8xl0 54 s _1 . The SED adopted here is 
somewhat different from the one used to model the X- 
ray absorption in Paper IV, which had a =0.5, 4.7, and 
0.77 in the energy ranges 0.002 - 0.04, 0.04 - 0.1, and 
0.1 - 50 keV). We recomputed our models derived in 
the next section with the SED used in Paper IV and 
find no significant differences in the predicted column 
densities for the ions measured in the UV spectrum, af- 
ter scaling the ionization parameter to account for the 
somewhat larger relative EUV flux used in the X-ray 
models. The issue of SED is addressed further below in 
the discussion of the combined UV and X-ray modeling 
results (§6.2 ). We adopted the rough ly solar elemental 
abundances l|Grevesse fc Anders 119891) used in the KC01 
models, and assumed no dust is present in the absorber. 

The combined UV and X-ray spectrum of NGC 3783 
shows that a large range of ionization states is present in 
all kinematic components and there are multiple zones of 
ionization overlapping at all absorption velocities (KC01; 
Paper IV). Thus some lines, particularly from more 
highly ionized species, may be comprised of blends of dif- 
ferent physical components having different covering fac- 
tors and optical depths; evidence for this was presented 
in §3 for kinematic component 1, for example. There- 
fore, we take as primary constraints the lines least likely 
to be affected by blending or saturation. These include 
lines from the lowest-ionization species detected in each 
kinematic component, ions with low abundances of the 
parent element (e.g., P V), and ions in excited states (e.g. 
C III*). Additionally, upper limits on column densities 
from non-detections provide unambiguous constraints. 

5.2. Model Results for Mean Low and High State 
Spectra 

Figure 7a shows the solutions in log(J7) and \og(Nu) 
from a grid of photoionization models that match the 
measured column densities for component 1. The con- 
tours on the left are for the low covering factor absorber, 
component la. The thickness of the contour for each ion 
spans the range of solutions corresponding to estimated 
uncertainties in t he measured column densities (e.g., see 
lArav et aLll200lD . Figure 7a includes solutions for Si IV, 
C II, C III, and P V measured in the low-state spectrum, 
plotted with hatched marks. The C III solutions are for 
the measured metastable level column densities, using 
the electron density derived from the C III* A1175 fea- 
ture in §4 and an electron temperature of T e wl.8xl0 4 
K, consistent with the best-fit Cloudy model. The solu- 



tions to the high-state column densities (mean of S15 and 
S16) for Si IV are also included on the plot, shown as the 
contour with no hatch marks. These were shifted along 
the x-axis by log(f/)= —0.23, corresponding to the peak 
amplitude of flux variation observed during the intensive 
monitoring. This gives the high-state solutions in terms 
of the low-state ionization parameters, allowing a direct 
comparison of the two flux states on the same grid: for a 
correct model with both states in ionization equilibrium 
and an ionization parameter that scales as the observed 
UV continuum flux, the low-state and shifted high-state 
solution contours would overlap. High-state solutions oc- 
cupying regions in U - Nh above and to the left of the 
low-state solutions overestimate the variability observed 
between states, while solutions to the right and below 
underestimate the variability. 

Figure 7a shows all ions plotted for the low-state spec- 
trum are fit well by an extended, narrow range of models 
in U - Nh parameter space, with lower bounds log{U) > 
— 1.7 and log(Nu) > 20.3. For higher ionization models, 
these ionic abundances (particularly Si IV and P V) are 
very sensitive to small changes in U and Nh, as seen 
in the narrowing of the solutions spanning the measured 
limits on the column densities. This is due to the He II 
opacity - these solutions are in the region in parame- 
ter space where the He II edge becomes optically thick. 
In these models, the relatively low-ionization species in 
component la exist primarily in a small region in the 
back end of the slab, where the continuum flux is heav- 
ily filtered, thus their abundances are greatly affected by 
the sensitivity of the He II opa city in the Stromgren shell 
to the model parameters (e.s. iKraemer et aTl t20021 . For 
these high-ionization models, the range of solutions be- 
comes linear in U - Nh, tracing the He II column density 
contour. 

Combining the shifted solution to the high-state Si IV 
column density, the range of solutions is limited to low- 
ionization values. The overlap between Si IV high and 
low-states gives slightly lower U than the fit to all low- 
state lines, but the solutions are close. For models with 
log(U) > —1.5 in the low-state, Figure 7a shows the 
Si IV variability is predicted to be significantly greater 
than observed. This discrepancy becomes pronounced for 
higher-ionization solutions: e.g., a factor of 1.7 increase 
in ionizing flux from the logf = —1, logiVff =21.6 model 
results in a Si IV column density of ~ 10 12 cm -2 , which 
is 40 times weaker than observed. Thus, based on these 
models of the mean high and low-state spectra, the lowest 
ionization absorber in component 1 has best-fit solution 
log(L0 = -1.6±£; 2 , log(iVff) = 20.6±°; 4 . The limits on 
the high-state C II, C III*, and P V column densities 
(not shown on Figure 7a) are also compatible with these 
solutions. We note the predicted column densities for 
O VI, N V, and C IV in this absorber would produce 
heavily saturated absorption in their resonance doublet 
lines, with Ni on > 5xl0 15 cm -2 for each ion. Model 
predictions for all ions are listed below their measured 
values in Table 2. 

The physical parameters of the absorber with higher 
covering factor, component lb, can also be constrained. 
For optically thin conditions, the C IV and N V column 
densities each scale linearly with Nh, but have different 
dependences on U; thus, the C IV : N V ratio uniquely 
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determines the ionization parameter, independent of the 
total column density. As a result, assuming the ioniza- 
tion structure in component lb is uniform over all radial 
velocities, the C IV and N V column densities measured 
in the unsaturated red wing determines U for the en- 
tire absorber. Our measurements give log(Z7)= — 0.4^Q2. 
A lower limit on the total column density comes from 
the lower limits measured on C IV and N V over the 
full, saturated profile. For the above U, this gives 
log(Vff) > 20.3. The O VII column density measured 
in the CXO spectrum over the velocity range coinciding 
with component 1 (see Figure 10 in Paper I) provides an 
upper limit on Nh ■ Our fit to this high- velocity region of 
the profile gives Novn = 5xl0 16 cm~ 2 , with an upper 
limit of < 10 18 cm -2 at the 90% confidence level. Incor- 
porating this as an upper limit on O VII gives Njj < 21 .4 
for component lb. The region in U, Nh parameter space 
spanned by these limits is also plotted on Figure 7a. 

The modeling constraints for component 2 are C IV 
and N V in both the low and high states, and the upper 
and lower limits on C III and H I, respectively. These 
solutions are shown on the U — Nh plot in Figure 7b; 
N V and C IV are shown as hatched contours in the 
low-state and as unfilled dashed contours in the high- 
state, which are shifted to account for the flux difference 
between states as in Figure 7a. In each state individ- 
ually, the C IV and N V solutions overlap over an ex- 
tended range in U — Nh', however, the combined high 
and low-states are simultaneously fit by only a small re- 
gion of parameter space selecting the lower ionization so- 
lutions; log(t/) = -0.45t°; 2 (low-state values), \og(N H ) 

= 20.41°'^. These solutions are also consistent with the 
H I and C III limits. 

In component 3, C III A977 is clearly detectable in 
the low-state FUSE spectrum, representing the lowest- 
ionization species detectable in this component, but has 
only an upper limit in the high state (see §3.2.5 and Fig- 
ure 4b). For our adopted effective covering factor, C IV is 
near saturation in the low-state, but well below the resid- 
ual flux level in the high-state giving a correspondingly 
small column density. The N V doublet is consistent with 
being saturated in both states. The model results are 
shown in Figure 7c, with high-state solutions shifted as 
in Figures 7a and 7b. This shows the C IV high-state and 
C III low-state column densities (hatched contours) are 
simultaneously matched by a restricted space of corre- 
lated U and N H values, with log(t/)> -1.2, log(AT ff )>19. 
Including lower limits on N V and H I places lower limits 
of \og(U) > —0.7 and > —0.6, respectively. An upper 
bound on the solution comes if we require the model X- 
ray columns do not exceed those measured in the CXO 
spectrum. The O VII column density integrated over 
all radial velocities {Nqvii = 10 18 cm~ 2 from Paper I), 
provides the most rigid constraint and is plotted in Fig- 
ure 7c. Incorporating this limit and assuming all the lines 
measured in the UV arise in a single physical component 
gives log(E0= r 0.5±g;i, \og{N H )= 21.1±g;l- 

These modeling results can be compared to those from 
KC01, which were based on the first STIS observation. 
The ionization and total column density for compo- 
nent la are substantially greater here than in KC01. This 
is due to the constraints implied by the detection of C II 
and P V in the data presented in this study, and the 



use of a lower covering factor for Si IV, giving a larger 
Si IV column density. A reliable covering factor for Si IV 
was made possible by the high S/N in the merged STIS 
spectrum (Paper II). The solutions to components 2 and 
3 presented above have similar total column densities as 
the KC01 solutions, but have somewhat lower ionization 
parameter (~ 0.3 - 0.4 dex). For component 2, this lower 
ionization was imposed on the solution primarily by the 
magnitude of variability observed in C IV and N V be- 
tween high and low flux states. For component 3, the 
detection of C III A977 in the FUSE spectrum and the 
limit placed by O VII measured in the CXO spectrum 
drove the solution to lower ionization. 

5.3. Time- Dependent Ionization Solutions 

Here we explore variability timescales for the UV ab- 
sorbers. With the density constraint derived for com- 
ponent la from C III*, the time-dependent ionic pop- 
ulations for this absorber can be probed in detail for 
comparison with the observed variability, thereby pro- 
viding better constraints on its physical conditions. For 
the other absorption components, the observed magni- 
tude and timescales for variability in ionic populations 
in response to the continuum flux variations can be used 
to place limits on the density. 

5.3.1. Detailed Variability Calculations for 
Component 1 

The response time for an ion is a strong function of 
several factors: its population (m) relative to adjacent 
ionization stages, the electron density, and the mag- 
nitude of change in ionizing fl ux incident on the ab- 
sorberfe.g. lHamann et all 119971) . We computed time- 
dependent ionic populations in response to changes in the 
ionizing continuum for the gas in component 1 based on 
the density and ionization solutions determined above. 
The problem involves solving for ionic abundances, rij, 
from the system of first-order differential equations: 

-rr = -[Pi+ n earec,i-i(T)}ni+72 e ni + ia rec ,i(T)+pi-ini-i, 

(6) 

where pi and n e a reCj i are the ionization rates from and 
recombination rates to ionization stage i (e.g., Krolik 
& Kriss 1995). Output from the low and high-state 
Cloudy equilibrium models were used for the param- 
eters in equation 6, corresponding to initial and final 
states associated with a change in ionizing flux. The 
time-dependent ionic abundances were then solved us- 
ing the Runge-Kutta-Fehlberg method, which compares 
fourth and fifth-order Runge-Kutta estimates to adjust 
the step size. The rates include all ionizati on and recom- 
binat ion processes treated by Cloudy fsee lFerland et af] 
1998). We consider simple step- function increases and 
decreases in the ionizing flux below. 

For example, to compute the time-dependent compo- 
nent 1 Si IV column density in response to an increase in 
ionization from the low to high states, the initial values of 
all ionic species of silicon (ni(t — 0)) were set to the val- 
ues from the low-state equilibrium model (\og(U)= —1.6, 
log(iVjj)=20.6). The ionization rates, pi, were set to 
the values from the corresponding high-state equilibrium 
model (log(?7)= —1.37). This assumes the ionization 
rates change instantaneously, which is valid if they scale 
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linearly with the ionizing flux, as is the case in the model 
considered here. The recombination rates differ some- 
what between the initial and final state equilibrium mod- 
els because of their temperature dependence. However, 
for the time intervals and magnitude of flux variations 
considered here, the absorber is not in thermal equilib- 
rium, as determined by the thermal timescale: 

_ 3/2n e kT 

Hhermal — n 2yWy) ' ^ > 

where the numerator is the total thermal energy and the 
denominator the net cooling rate in the gas per unit vol- 
ume, with A(T) representing the difference between to- 
tal heating and cooling rates, which are obtained from 
the Cloudy calculations. Even for extreme changes in 
ionizing flux, the thermal timescale (taken here as the 
e-folding time) is of order 1/2 year for these conditions. 
It is much longer than this for the relatively minor flux 
changes observed in NGC 3783, due to the small value 
of the net cooling rate. Thus, the gas is not in ther- 
mal equilibrium and the actual value of the temperature 
reflects an average of the long-term history of the cool- 
ing and heating rates as they respond to variations in 
the ionizing flux. For the conditions considered here, 
the temperature can be assumed to be constant. Thus, 
we adopt the recombination rates from the initial state 
models in our solutions to equations 6. 

Figure 8 illustrates the dependence of time-dependent 
populations on the amplitude of ionizing flux variations. 
The time-dependent behavior of the Si IV column den- 
sity is given for both increases and decreases in ioniz- 
ing flux by factors 1.4, 1.7, and 3, showing the time 
needed to achieve a given level of variability is sensitive 
to the amplitude of flux change. For models with in- 
creased ionization, we used the low-state model solution 
as the initial state, and for decreased ionizing flux, the 
high-state model. For large decreases in ionizing flux, 
the recombination timescale (i.e. , e-folding time) is see n 
to approach the expression from lKrolik fc Krissl Jl995), 
trecomb = rii / (n i+1 n e a rec ^); t recomb = 5 days for Si IV 
for the model considered here. Figure 8 also shows the 
time interval corresponding to a given fractional change 
in the ionic abundance is similar for cases of ionization 
and recombination for the moderate flux variations con- 
sidered here. However, for very large flux changes, our 
calculations show the ionization timescale becomes much 
shorter than the recombination timescale. 

We calculated the detailed variability for the compo- 
nent la absorber based on the observed flux variations in 
individual STIS observations. In Figure 9, calculations 
for multiple step-function variations in flux, correspond- 
ing to the continuum light curve during our intensive 
monitoring, are compared with the measured Si IV col- 
umn densities. The solid line shows results for the model 
derived in §5.2 based on the low and high-state solutions; 
the UV light curve is shown in the bottom panel for com- 
parison. The rates and initial ionic populations were 
taken from the low-state Cloudy model, log(i7)=— 1.6, 
logiVff =20.6. At each epoch, the ionization rates were 
then scaled by the change in flux observed in the subse- 
quent observation, and the time-dependent populations 
computed using equation 6. The overall Si IV variability 
is seen to be matched well by this model. It reproduces 
the somewhat damped decrease in Si IV column density 



observed in high-states S10, S12 following the increased 
flux in those observations, and the further decrease in 
later epochs, where the absorption is a minimum at S15. 

In §5.2, an extended range of solutions, with correlated 
U and Nh values, was found to match the low-state col- 
umn densities in component la reasonably well, with the 
high-state solution selecting the lower U models (Figure 
7a) . If the physical conditions in this absorber are such 
that the Si IV population is not able to respond to the 
ionizing continuum sufficiently, the gas could in principle 
be in a higher ionization state than modeled by assum- 
ing full equilibration between states. With the density 
for component la determined independently from C III*, 
it is possible to test this directly. Thus, we computed 
the time-dependent Si IV population for higher U mod- 
els to test its variability. Figure 9 shows results for a 
model with \og(U) = -1, log(N H ) =21.5 (dashed line), 
which was selected from the low-state solutions in Fig- 
ure 7a. This high U model is seen to far overestimate the 
variability in Si IV, and thus is excluded by our timing 
analysis. 

5.3.2. Constraints on Density in Components 2 and 3 
from Observed Variability 

Components 2 and 3 do not have direct estimates of 
the density, and thus are less well constrained than com- 
ponent 1. However, lower limits on their densities can 
be derived based on the observed variability. This con- 
straint comes from the characteristic time separation be- 
tween the low and high-state epochs, which provides an 
upper limit on the timescale for changes in ionic abun- 
dances as they respond to the ionizing continuum. 

For component 2, we computed the time-dependent 
C IV and N V abundances in the same manner as above 
for component 1 to derive these limits. We assumed a 
simple step- function change in flux of a factor 1.7, corre- 
sponding to the continuum variation between the mean 
low and high-states during the intensive monitoring. We 
parameterized the observed variability for these mod- 
els by requiring the ionic populations varied by at least 
half the amount observed over a characteristic timescale, 
which we define as the time interval between the mean of 
the low states (S5 - S8) and mean of the high states (S10, 
S12), ~ 20 days, where maximum change in these column 
densities occurred. This places a lower limit on n e . The 
c-folding variability times are proportional to the inverse 
of the density (t oc n~ ). Thus, we computed the model 
for a fiducial density and then solved for the value of n e 
that reproduces the required level of variability over the 
characteristic timescale. 

To demonstrate this, Figure 10 shows the time- 
dependent N V and C IV ionic abundances for com- 
ponent 2 based on the solution in §5.2 (log([/)=— 0.45, 
log(A r n)=20.4). These were computed from equation 6 
for a fiducial density of n e =10 4 cm~ 3 , thereby giving 
Ni on /Ni on j=o as a function of t x (n e /10 4 ). Horizontal 
dashed lines mark the 50% level of variability; dotted 
lines give the final, equilibrium values. The lower lim- 
its on density follow straightforwardly from the intersec- 
tion of the abundance curves with the variability limits 
(dashed lines), which is assumed to occur over a 20 day 
interval: n e — t/20x 10 4 cm~ 3 . The N V variability gives 
n e > 1.3xl0 3 cm" 3 and C IV gives n e >2.3xl0 3 cm" 3 
for component 2. Component 3 is less well constrained 
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because for each of the line constraints, at least one state 
provides only a limit on the column density. A similar 
analysis for this component gives n e >7.5xl0 2 cm~ 3 . 

6. INTERPRETATION 

6.1. Constraints on Physical Conditions and Geometry 
of the UV Absorbers 

We now derive constraints on the physical conditions 
and geometry of the absorbers based on our above anal- 
ysis. The best constraints are for the high-velocity out- 
flow region, component 1, due to the density measured 
directly from the metastable C III* feature (§4). From 
the expression for the ionization parameter given in §5.1, 
the distance between the absorber and central source 
can be solved for the low-ionization absorber (compo- 
nent la), using the derived values of the number density 
{nn ~ n e /1.2), Q, and U from the modeling in §5.2. 
This gives R =7.7^ x 10 19 cm (25 pc). If the individ- 
ual kinematic components are uniform clouds, i.e., hav- 
ing internal volume filling factors of unity, then their ra- 
dial physical depths can be determined straightforwardly 
from the ratios of the total column density to the num- 
ber density, giving AR =1.3xl0 16 cm for component la. 
Results are summarized in Table 3. 

In Paper III, we showed the radial velocities for all com- 
ponent 1 lines in the STIS spectra decreased at a rate of 
~ 50 km s _1 yr _1 . This includes Si IV, which comes from 
the low-ionization, low-covering subcomponent, and the 
N V and C IV lines which are predominately from com- 
ponent lb. Thus, these absorbers appear to be dynam- 
ically linked and, hence, co-located. Using this result, 
and adopting the ionization parameter for component lb 
derived in the unsaturated red-wing of C IV and N V 
(§5.2), the density for this higher ionization component 
can be solved from: nn,ib = n-H,ia x Uia/Uib, giving 
nn,ib = 1.9xl0 3 cm -3 . A broad range of values for the 
radial depth for component lb follow from the loose con- 
straints on N H , 8xl0 16 < AR <lxl0 18 cm. This can be 
compared to the lower limit on the projected transverse 
size of this absorber, determined by the partial coverage 
of the BLR from the Lyman line analysis and the size 
of the BLR based on reverberation mapping (Onken & 
Peterson 2002), giving X T > lxlO 16 cm (Paper II). The 
transverse size of component la is not well constrained 
because the details of the covering factors of the individ- 
ual sources are not known; the only constraint is that 
the UV continuum source is at least partially covered by 
this absorber. 

For components 2 and 3, the constraints derived on 
n e from the observed variability (§5.3.2) give upper lim- 
its on their distances from the AGN and radial depths. 
Our derived lower-limit on n e in component 2 based on 
the calculated response times for C IV and N V gives 
R < 7.4xl0 19 cm (24 pc) and AR < 10 17 cm. The 
density limit for component 3 implies R < 1.4 xlO 20 cm 
(45 pc) and AR <1.3xl0 18 cm. Lower limits on the 
projected transverse dimensions of these absorbers fol- 
low from the BLR covering factors, giving Xt > 1 - 
2xl0 16 cm. Noting that AR oc njj 1 and R oc njj 1 ^ 2 , 
the absorber's geometries are seen to depend strongly 
on their locations when combined with the constraints 
on the transverse sizes. For example, if the absorbers 
are uniform clouds having roughly spherical geometries, 



with dimensions approximately equal to the lower limits 
given by the BLR covering factors, then distances ~ 10 
and 5 pc are implied for components 2 and 3, respec- 
tively. Smaller distances would require flattened geome- 
tries. For example, if they are located at 1/2 pc, then 
they would be extremely thin shell structures, having ra- 
dial dimensions at least 100 - 500 times smaller than 
their transverse sizes. 

The results summarized in Table 3 show all UV kine- 
matic components have low global volume filling factors 
(AR/R <C 1), and are consistent with being relatively 
small, discrete clumps. Additionally their derived dis- 
tances (or limits) are consistent with the location of the 
inner narrow line region in AGNs, and inside the more 
extended, diffuse NLR. The geometry implied by these 
results is consistent with our earlier assumption that the 
NLR is unocculted by the individual components of UV 
absorption in deriving the covering factor model (§3.2.2). 

6.2. Connection between the UV and X-ray Absorption 

Here, we explore the connection between the observed 
UV and X-ray absorption. The absorption in the two 
bandpasses shows similar kinematic structure: in Papers 
I and II, all X-ray lines having sufficiently high-resolution 
and S/N were found to span the radial velocities of the 
three UV kinematic components. Modeling of the CXO 
spectrum in Paper IV showed the X-ray absorption is 
highly inhomogeneous, requiring three ionization zones 
that span a range > 50 in U to reproduce the full set of 
lines. A summary of the modeling and geometric con- 
straints from that analysis is given in Table 3, together 
with our results for the UV absorbers, which imply fur- 
ther inhomogeneities in the outflow in NGC 3783. To 
convert the X-ray component solutions for a consistent 
comparison with the UV models, we matched models de- 
rived with the SED defined in §5.1 with the model pre- 
dictions in Paper IV, which used an SED with stronger 
relative emission in the EUV (see §5.1). The ionization 
parameters were scaled by matching the peak ions in each 
ionization component with the predicted columns in Ta- 
ble 3 of Paper IV. The models computed with the two 
SEDs were found to give very similar results, matching 
the dominant ions in each component to better than 90% 

Table 3 shows the three UV components with relatively 
high-ionization (lb, 2, and 3) share the same ionization 
parameter as the lowest-ionization region modeled in the 
X-ray (hereafter XLI), although they have a smaller (in- 
tegrated) total column density. In Table 4, the predicted 
column densities from the UV absorbers for key ions with 
lines in XLI are listed. This shows the component 3 and 
component lb absorbers may give significant contribu- 
tion to some of the XLI lines. Indeed, upper bounds on 
U, Nh in these components are from the measured O VII 
column density (§5.2). Table 4 shows component 3 may 
also contribute strongly to Mg IX and about a third of 
the Si IX and Si X measured in the X-ray spectrum. If 
Nh is near the upper limit for the less well constrained 
component lb absorber, it may produce strong Si IX - 
Si XI, Mg VIII, and Mg IX. O VI is the only ion with 
lines detected in both bandpasses. The UV doublet is 
found to be saturated in all three kinematic components 
(Figure 4b), giving only a lower limit, and thus con- 
sistent with the measurement from the CXO spectrum, 
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Novi >10 17 cm~ 2 (Paper IV). 

The XLI model predicts a C IV column density that 
would produce heavily saturated absorption in the UV 
doublet in all flux states observed during our monitor- 
ing. However, as shown above, the C IV equivalent width 
varied in components 2 and 3, implying these features 
are not saturated, at least in the high flux states when 
their absorption was weaker. Although C IV is likely 
saturated in component 1, the dominant X-ray absorp- 
tion coincides kinematically with the lower velocity UV 
components (Paper I). Our model for component 3, con- 
strained by the unsaturated high-state C IV, fails to re- 
produce the lowest ionization X-ray species measured in 
Paper IV, underestimating Si VII, Si VIII, and O V by 
factors of 10 or more. Component 2 contributes even less 
due to its relatively small total column density. We have 
done extensive modeling to test if this discrepancy be- 
tween the measured UV and X-ray column densities can 
be reconciled for any values of the model parameters. We 
find no reasonable choice for an SED that reproduces the 
measured Si VII and O V column densities while main- 
taining unsaturated absorption in the C IV UV doublet, 
due to the similarity in ionization potentials of these ions, 
nor does any combination of multiple components with 
different physical conditions. One possible explanation 
for this discrepancy is that it is due to a geometrical 
effect. For example, the very large column density of 
C IV implied by the XLI model may be buried in the 
variable, unsaturated absorption in components 2 and/or 
3. This would require a low covering factor of the UV 
emission by XLI given the relatively shallow absorption 
observed in C IV components 2 and 3 (see Figure 4a). 
Since the X-ray emission region is much more compact 
than the UV BLR (and likely the UV continuum source 
as well), it is reasonable that the same gas would have 
different line-of-sight covering factors in the two band- 
passes. Qualitatively, this is consistent with the inho- 
mogeneous wind model described in detail in the next 
section, in which denser, lower-ionization regions occupy 
smaller and smaller volumes in the global outflow. The 
lowest ionization X-ray species could conceivably come 
from a region that is sufficiently dense and compact that 
its presence is not seen against a larger absorber giving 
the unsaturated C IV. Alternatively, it may imply the 
XLI absorber does not occult the UV absorber at all and 
that we are seeing physically distinct absorbers in dif- 
ferent lines-of-sight in the two bandpasses, thus placing 
very specific requirements on the absorption - emission 
geometry. 

Finally, we compare independent constraints on the 
distances of the absorbers based on the UV and X-ray 
analysis. In Paper IV, limits on variability in the X-ray 
absorption were used to give lower limits of R > 3, 0.6, 
0.1 pc for XLI, XMI, and XHI, respectively (see Table 
3). Similar lower limits were determined from variability 
analysis of XMM-Newton observations by Behar et al. 
(2003) . In contrast, a recent study bv iKrongold et all 
(2005) reported variations in the Fe M-shell unresolved 
transition array in the CXO observations, and they de- 
rive an upper limit of « 6 pc for this absorber. Limits 
can also be derived on the distances for each of the ion- 
ization components modeled in Paper IV based on the 
simple geometrical requirement that AR < R. If the ab- 
sorbers are uniform, constant density regions, then for a 



given model solution with parameters U - Nh , AR oc R 2 
through their dependences on density. This gives up- 
per limits on R (lower limit on tih); results are listed 
in Table 3. Filling factors <1 would imply more strin- 
gent limits on R, since the absorbers would then occupy 
a larger region than given by AR. The limit on XHI, 
R < 4 pc, is somewhat less than the distance derived 
for component la based on the C III* density constraint 
and UV modeling above; it is similar to the estimates for 
UV components 2 and 3 based on the assumption of uni- 
form absorbers, with roughly spherical geometries (see 
discussion in §6.1). 

6.3. Global Model of the Outflow m NGC 3783 
6.3.1. Evidence for an Inhomogeneous Wind 

With constraints derived on the physical state and ge- 
ometry of the absorbers, we now investigate implications 
for the global model of the outflow in NGC 3783 and ex- 
plore constraints for dynamical models. In Paper IV, the 
three ionization components modeled for the X-ray ab- 
sorption were found to be consistent with being in pres- 
sure equilibrium (P oc T/U), all lying on stable regions 
of the nearly vertical part of the thermal stability curve 
(log(T) vs log(E//T); see Figure 12 in Paper IV). Based 
on the modeling above, UV kinematic components lb, 2, 
and 3 are also at the same pressure. They occupy the 
low-temperature base of the region of the thermal sta- 
bility curve where a range of temperatures can co-exist 
in pressure equilibrium, coinciding with the solution in 
log(T), log(C//T) for the low ionization X-ray gas, XLI, 
seen in Figure 12 in Paper IV. Thus, these results are con- 
sistent with the general model presented for NGC 3783 in 
Paper IV, an d described theoretically in lKrolik fc K riss 
1)19951 1200 lj) . in which the absorption arises in a multi- 
phase thermal wind, comprised of embedded regions that 
are inhomogeneous in temperature and density. In this 
model, the UV absorbers (and XLI) represent the lowest- 
ionization, densest material detectable in the spectrum. 

These results can be compared furt her with the recent 
study bv IChclouch e fc Netzer I l)2005fl . which presents a 
new dynamical model for the X-ray outflow in NGC 
3783 that combines some of these general ideas with 
more specific assumptions and calculations. The main 
driver of the flow in this model is thermal gas expan- 
sion and the main carrier of the flow is the highest ion- 
ization, hottest component. The model suggests that 
cooler, lower-ionization material occupies smaller frac- 
tions of the flow where the size distribution resembles 
what is known from the ISM. An inhomogeneous wind 
model opens the possibility for having different covering 
fractions for different ionization components that share 
the same outflow velocity. It also eases the problem 
of the large transverse dimension of the X-ray gas in 
component XLI discussed in Paper IV, since the low- 
ionization components can have large column densities 
yet they are composed of small filaments or clouds with 
relativ ely small dimensions. In th e framework of this 
model, iChelouche fc Netzer I (|2005l) compute the kinetic 
energy associated with the NGC 3783 outflow to be only 
a small fraction of the bolometric luminosity, and the 
mass-loss rate to be comparable to the accretion rate. 

The low-ionization, high-velocity UV absorber compo- 
nent la does not fit into the picture of inhomogeneities 
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co-existing at pressure equilibrium. It has a gas pressure 
that is a factor of ten greater than the other compo- 
nents and thus, if embedded in the more diffuse higher- 
ionization gas without an additional confining mecha- 
nism, will eventually evaporate. One possibility is that 
component la is comprised of relatively high density ma- 
terial that has recently been swept up from an external 
mass source and exposed to the ionizing radiation from 
the AGN, and is destined to expand to come into pres- 
sure equilibrium with the remainder of the flow. This 
absorber was found to have the smallest covering fac- 
tor, based on the Si IV doublet absorption, consistent 
with this being a dense, compact region in the flow. Per- 
haps it is embedded in, and evaporating into the more 
diffuse component lb gas, which is at pressure equilib- 
rium with the rest of the outflow. This would explain 
the dynamical link between these regions implied by the 
decrease in radial velocity observed in both absorbers 
(Paper III). If component la is unconfined, the timescale 
for it to expand, with its density decreasing to that of 
the lb region, can be approximated based on the es- 
timated size of the absorber (AR « 2xl0 16 cm) and 
the thermal velocity (vth = 22 km s _1 at the model 
temperature of 1.9xl0 3 K), giving a lifetime of ~ 150 
years. Alternatively, component la (and all UV ab- 
sorbers) could be confined by magnetic pressure, as pro- 
posed in some dynamical model s (e.g. lEmmering et alJ 
Il992tlde Kool fc Beeelman Il995|) . Only a moderate field 
(B ~ 10~ 3 G) would be required to balance the t hermal 
pressu re of this absorber; general calculations bv iRees 1 
( 1987) show fields of this strength could easily be present 
at the distance derived for the UV absorbers in NGC 
3783. 

The independent appearance of the UV components on 
yearly timescales, without correlation with the observed 
continuum flux (KC01), provides further constraints for 
physical models of the outflow. In the framework of the 
thermal wind model described above, it may be a signa- 
ture of condensations in the outflow that have appeared 
in our linc-of-sight to the AGN. However, we note these 
timescales for the appearance of absorption components 
are quite short compared with the evaporation timescale 
derived for component la above. Alternatively, it may be 
due to motion of the UV absorbers across our sightline, 
as discussed in KC01 and Paper III. In this scenario, dy- 
namical models must account for the implied transverse 
component of velocity (v > 500 km s _1 , KC01; Paper 
III), at the ~ 10 pc distance scale, as well as the radial 
velocity for the UV absorbers. Additionally, the decreas- 
ing radial velocity observed in component 1 must be ex- 
plained, consistently for both subcomponents la and lb. 
If this is due to a geometrical effect, in which the ab- 
sorber is following a curved path across our line-of-sight, 
it provides constraints on its trajectory and kinematics 
(see Paper III). 

Finally, the combined UV and X-ray results on the 
outflow geometry and kinematics may provide additional 
constraints. If all ionization components are roughly co- 
located, then the gas giving rise to the UV absorbers 
would occupy a much smaller region than the high- 
ionization material seen in the X-ray. If the UV ab- 
sorption components are single uniform clouds, it may 
be difficult to account for the similarity in outflow ve- 
locities for the full range of ionization seen in the two 



bandpasses; it would require three small, localized re- 
gions (the UV absorbers) span the full kinematic range 
of the much more extended global outflow correspond- 
ing to the higher ionization gas. One possibility is that 
the individual components have low internal filling fac- 
tors, so that the lower-ionization gas we see is distributed 
over a larger volume, rather than localized to the small 
region implied by AR derived above. Another potential 
explanation is that the gas is not co- located (e.g., see dis- 
cussion in §6.2). This could be explained, for example, 
by dynamical models that predict the terminal veloci- 
ties of flow regions will depend on the launch radius of 
th e outflow, such as th e MHD wind models presented 
in iBottorff et al.l ((2000); in this framework, it is possi- 
ble that gas with a range of physical conditions and at 
different locations would have similar terminal velocities. 

6.3.2. Connection to the Observed Emission- Line 
Spectrum 

Here we explore a potential connection between the 
UV absorbers and the optical - UV line e mission. In 
a recent study, iCrenshaw fc Kraemer I l)2005|) argued the 
intrinsic UV absorbers in Seyfert galaxies may be iden- 
tified with a high-ionization component of gas in the 
inner narrow line region (NLR), based on their kine- 
matics, ionization, inferred distance, and global cover- 
ing factor. This is supported by our results for NGC 
3783. Fr om STIS slitless spectra of a sample of Seyfert 
galaxies. ICrenshaw fc Kraemer I (|2005|) found 9 of 10 ob- 
jects have bright, central compact [O III] emission-line 
knots with half-width at half-maximum sizes of ~ 5- 
70 pc and half-width at zero intensity radial velocities 
of 300-1100 km s _1 . The derived distances and ob- 
served outflow velocities for the intrinsic UV absorbers 
in NGC 3783 are consistent with this inner NLR com- 
ponent. A further comparison can be made with the 
line-emission obse rved in NGC 3 783. Comparing to the 
measurements in IE vans I (^988), we find the UV ab- 
sorbers may contribute significantly to some of the high 
ionization UV-optical emission-line fluxes. For compo- 
nent la, the most constraining line is [Ne V] A3426, 
with a luminosity of L salO 41 e rgs s~\ corrected for th e 
NLR reddening determined bv IWard fc Morris I (l984). 
This can be compared to the prediction from compo- 
nent la, giving L =3xl0 41 ergs s _1 xC g , where C g is 
the global covering factor of the absorber. This limits 
C g < 0.3 for gas with similar physical conditions as com- 
ponent la. This component also predicts a strong contri- 
bution to the [O III] A5007 emission, and the coronal iron 
lines ([Fe VI] A5177, [Fe VII] A6085, [Fe X] A6375, and 
[Fe XI] A7891), reproducing 30 - 100% of their observed 
line luminosities for C g =1. Absorption component 3, 
with its relatively high-ionization and large column den- 
sity, is also predicted to contribute strongly to the highcr- 
ionization coronal iron lines; e.g., it produces ~3 times 
the observed [Fe XI] A7891 luminosity for full global cov- 
erage, thereby limiting C g <0.3. It may also contribute 
to the observed high-ionization UV resonance lines, par- 
ticularly O VI AA1032,1037, implying a similar limit on 
C g , with only a small contribution (w 5 - 10%) coming 
from resonance scattering. These constraints on C g for 
the absorbers in NGC 3783 are similar to the general 
constraint implied by th e detection rate of UV absorp- 
tion in Seyfert Is f~50%. ICrenshaw et al.lll999|) . Finally, 
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we note there is potential evidence for the connection be- 
tween the UV absorbers and high-ioniza tion line emission 
in the profiles observed in NGC 3783. iWard fc Morris I 
( 1984) found a strong asymmetry towards blue wave- 
lengths in the forbidden coronal iron lines, with emission 
extending to w —1500 km s _1 , which is not present in the 
low-ionization NLR lines such as [S II] (see their Figure 
6). 

7. SUMMARY 

We have presented analysis of the intrinsic UV ab- 
sorption in the Seyfert 1 galaxy NGC 3783, based on 
an intensive monitoring campaign with HST/STIS and 
FUSE. Our 18 STIS observations included both an inten- 
sive phase (3-8 day sampling) and a long-term phase 
(several months sampling), observing the UV continuum 
in a range of flux states spanning a factor ~ 2.5. In the 
three kinematic components with strong UV absorption 
(components 1, 2, and 3 at radial velocities v r ~ —1350, 
—550, and —725 km s _1 , respectively), the equivalent 
widths of the lowest-ionization lines present were ob- 
served to vary inversely with the UV continuum flux. 
This indicates all components are responding to the ion- 
izing flux, including the short timescales of the intensive 
observations. 

We isolated the NLR emission using variability in the 
emission-line profiles, in order to explore its effect on 
the line-of-sight absorption covering factors. We find a 
coverage model with an unocculted NLR and individual 
covering factors of the continuum and BLR sources de- 
rived from the Lyman lines predicts saturation of several 
lines that have substantial residual fluxes. These results 
are consistent with the independent evidence for satu- 
ration in these lines from the lack of observed variabil- 
ity and predictions from modeling results. This coverage 
model is used to obtain ionic column densities in low and 
high state spectra for photoionization modeling. Based 
on differences in derived covering factors and kinematic 
structure, the absorption associated with component 1 
appears to be comprised of two physical components (la 
and lb). The decrease in radial velocity observed in lines 
associated with both regions implies they are linked dy- 
namically and co-located. 

We were able to resolve the C III* A1175 absorption 
complex observed in component la into its individual 
multiplct lines and measure column densities for each 
metastable level. Comparing the ratios of the level pop- 
ulations with our calculations, we determine the electron 
density for this absorber to be ~2.5xl0 4 cm~ 3 . 

We compare the measured ionic column densities with 
photoionization models to constrain the physical condi- 
tions in each UV component, incorporating variability 
observed between low and high states. Based on the 
density from C III* and the modeling solution, compo- 
nent la is found to be located at ~ 25 pc. It is in a 
relatively low-ionization state (\og(U)= —1.6) compared 
with the other UV components. Using the density de- 
rived for this absorber, we compute time-dependent ionic 
populations and are able to reproduce the detailed vari- 
ability observed in this component. Components lb, 2, 
and 3 are found to have similar ionization parameters 
(log(Z7)= —0.5 0.4), with total column densities span- 
ning a range of a factor ~5. The ionization parameter 
for these components is the same as that modeled for 



the lowest-ionization X-ray region from the CXO spec- 
trum, but with smaller total column density. This may 
indicate different covering factors in the two bandpasses, 
perhaps due to the different sizes and line-of-sight geom- 
etry of the background emission regions. Analysis of the 
observed variability in components 2 and 3 limits their 
distances from the central ionizing source to be < 25 and 
< 50 pc, respectively. 

The modeling results for the higher ionization UV com- 
ponents (lb, 2, 3) imply pressures similar to those of the 
three X-ray ionization components modeled in Paper IV. 
Thus, the global outflow in NGC 3783 is consistent with 
an inhomogeneous wind, with relatively small dense re- 
gions giving the UV and lower-ionization X-ray absorp- 
tion, which are embedded in a hotter, more tenuous flow. 
However, component la has a factor of 10 greater pres- 
sure than the other absorbers, thus, if embedded in a 
larger outflow and unconfined it will eventually evapo- 
rate, over a timescale of ^150 years. Alternatively, this 
absorber could be confined by another mechanism, such 
as magnetic pressure. The emission-line luminosities pre- 
dicted from our models of the UV absorbers indicate they 
may be contributing substantially to the high-ionization 
lines observed in the UV-optical spectrum. The implied 
constraints on the global covering factors of the UV ab- 
sorbers in NGC 3783 are similar to those implied by the 
detection rate of intrinsic absorption in Seyfert 1 galax- 
ies. Additionally, the distances derived for the UV ab- 
sorbers are consistent with the location of the inner nar- 
row line region. 
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TABLE 1 
Archival UV Spectra of NGC 3783 



Instrument/ 


Coverage 


Resolution 




Date 




Grating 


(A) 




(A/AA) 




(UT) 




IUE SWP 


1150 


- 1978 




~250 




1978 November 26 - 


1992 July 29 a 


FOS G130H 


1150 


- 1605 




~1200 




1992 July 27 




GHRS G160M 


1527 


- 1597 




20,000 




1993 February 5 




GHRS G160M 


1527 


- 1597 




20,000 




1994 January 16 




GHRS G160M 


1527 


- 1597 




20,000 




1995 April 11 




STIS G140L 


1150 


- 1715 




~1000 




2000 May 17 




a See the IUE Merged Log at http://archive.stsci.edu/iue 


















TABLE 2 










Measured and Predicted Ionic Column Densities in the UV Absorbers'* 




Ion 


Component 1 




Component 2 


Component 3 




Low 


High 


Low 




High 


Low 


High 


H I 


>350 


>240 


>20 




>15 


>80 


>55 




(670) 


(420) 


(24) 




(12) 


(120) 


(66) 


N V 


>11 


>11 


11 ±3 




4.9 ±1.5 


>17 


>n 




(160) 


(160) 


(9.5) 




(3.5) 


(60) 


(25) 


N III 


>8 


4 


<0.5 




<0.5 


<1.5 


<1.0 




(34) 


(10) 


(-) 




(-) 


(0.05) 


(0.01) 


C IV 


>4.5 


>4.5 


2.2 ±0.5 


0.82 ±0.2 


>10 


4.1 ±1.4 




(340) 


(190) 


(2.7) 




(0.71) 


(21) 


(6.1) 


C III 






<0.5 




<0.5 


0.60 ±0.20 


<0.3 


C III* b 


(82) 


(25) 


(-) 




(-) 


(0.38) 


(0.08) 


0.8 (±0.25) 


0.35 (±0.3) 


<0.2 




<0.2 


<0.2 


<0.2 




(0.3) 


(0.08) 


(-) 




(-) 


(-) 


(-) 


C II 


0.48 ±0.2 


<0.3 


<0.25 




<0.25 


<0.10 


<0.10 




(0.59) 


(0.12) 


(-) 




(-) 


(-) 


(-) 


O VI 


>14 


>14 


>12 




>12 


>51 


>51 




(110) 


(51) 


(2.2) 




(1.0) 


(1300) 


(640) 


Si IV 


2.0 ±0.5 


0.55 ±0.25 


<0.1 




<0.1 


<0.05 


<0.05 




(1.2) 


(0.23) 


(-) 




(-) 


(-) 


(-) 


P V 


0.13 ±0.04 


0.09 ±0.06 


<0.2 




<0.2 


<0.2 


<0.2 




(0.12) 


(0.06) 


(-) 




(-) 


(-) 


(-) 



a Column densities listed are in units 10 14 cm 2 . All values measured using the 3-C covering factor model, except lines associated with 
component la (Si IV, C III*, C II, and P V), see text. Predicted values from best-fit models are listed below each measured column 
density. • ■ • denotes no measurement available due to blending with other features, ( - ) listed for model values <10 12 cm -2 . Lower 
limits are given for saturated lines. 

b Population of C+ 2 3 P° J =2 metastable level. 



TABLE 3 

Constraints on Physical and Geometrical Parameters in UV - X-ray Absorbers 



Component 


Log(C7) 


Log(7V H ) 


hog(n H ) 


Log(.R) 


Log(Afi) a 


Log(V T ) 






(cm~ 2 ) 


(cm -3 ) 


(cm) 


(cm) 


(cm) 


UV Component la 


-1.6 


20.6 


4.40 


19.94 


16.2 




UV Component lb 


-0.4 


20.3 - 21.4 


3.20 


19.94 b 


17.1 - 18.2 


>16.0 


UV Component 2 


-0.45 


20.4 


>3.3 


<19.9 


<17.0 


>16.3 


UV Component 3 


-0.5 


21.1 


>2.9 


<20.2 


<18.2 


>16.2 


X-ray LI C 


-0.5 


21.9 


0.6 - 4.7 


19.0 - 21.3 


>17.2 




X-ray MP 


0.8 


22.0 


2.0 - 5.0 


18.3 - 19.9 


>17.0 




X-ray HP 


1.3 


22.3 


3.2 - 5.4 


17.7 - 19.1 


>16.9 ■■■ 





a Radial depths assuming uniform, constant density absorbers. 

b Assumed co-located with component la based on decrease in radial velocities of lines in both components, see text. 

c Results for the three ionization X-ray components from modeling the CXO spectrum in Paper IV; U values converted for models with 
the SED used in the present study. Upper (lower) limits on njj (R, AR) based on limits on variability from analysis in Paper IV; lower 
(upper) limits based on AR < R. 
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TABLE 4 

Predicted Ionic Column Densities for Lines in X-ray SPECTRUM a 



Ion 


Comp la 


Comp lb b 


Comp 2 


Comp 3 


Total UV C 


cxo 


O V 


17.0 


15.1 - 16.3 


15.4 


16.2 


17.1 


17.1 (±0.40) 


O VI 


16.8 


15.9 - 17.2 


16.2 


17.0 


17.3 - 17.5 


>17.0 


O VII 


16.4 


16.8 - 18.1 


17.1 


17.8 


17.9 - 18.3 


18.04 (±0.20) 


O VIII 


14.5 


16.4 - 17.6 


16.7 


17.2 


17.4 - 17.8 


18.63 (±0.25) 


Si VII 


15.7 


13.4 - 14.7 


13.7 


14.7 


15.8 


16.26 (±0.20) 


Si VIII 


15.4 


14.5 - 15.7 


14.8 


15.6 


15.9 - 16.1 


16.66 (±0.08) 


Si IX 


14.6 


15.1 - 16.3 


15.4 


16.1 


16.2 - 16.5 


16.66 (±0.06) 


Si X 


13.5 


15.3 - 16.5 


15.6 


16.1 


16.3 - 16.7 


16.83 (±0.05) 


Si XI 


12.1 


15.1 - 16.3 


15.4 


15.8 


16.0 - 16.5 


16.29 (±0.06) 


Mg VIII 


15.1 


14.9 - 16.1 


15.2 


16.0 


16.1 - 16.4 


16.58 (±0.05) 


Mg IX 


14.1 


15.3 - 16.4 


15.6 


16.3 


16.4 - 16.7 


16.39 (±0.06) 



a Component 1-3 models are intermediate states between the low and high states used in modeling each UV component; all values are 
log(N lon ). 

b Range in ionic column densities corresponding to range in Njj for component lb. 

c Sum of model predictions from all UV components. 

d Column densities measured in the CXO spectrum, from Paper IV. 
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Fig. 1. — (la) Continuum fluxes measured at 1470 A spanning ~22 years of UV observations. Shown are observations from IUE, GHRS 
and FOS (A), and STIS (x) with error bars. The intensive IUE (JD » 2,448,500) and STIS (JD m 2,452,000) observations show variations 
of a factor of ~1.5 - 2.5 are typical on weekly/monthly timescales. The STIS observations span a range of moderately-low — moderately-high 
flux states, (lb) Continuum fluxes for observations presented in this study. STIS observations (x) are labeled chronologically SI— S18; 
fluxes extrapolated from the FUSE observations are plotted with (+) for the four FUSE observations that arc non-simultaneous with STIS 
observations (Fl, F3, F5, F6). The error bars represent estimated la uncertainties. 
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Fig. 2. — Intrinsic absorption equivalent widths as a function of UV continuum flux from all STIS observations. (2a) In component 
1, Si IV exhibits a general trend of decreasing absorption strength with increasing continuum flux, while C IV and N V show no strong 
variability. (2b) C IV and N V in component 2 are both variable, showing an inverse correlation with continuum flux. (2c) The C IV 
absorption in component 3 shows the same inverse correlation with continuum flux, whereas N V does not vary significantly. Observations 
during the intensive monitoring phase arc plotted with (o) to distinguish them from the long-term monitoring (x), and all observations 
are labeled chronologically for comparison with Figure lb. Results of least squares fits to all epochs are plotted for each line. The ratio of 
the computed slope to the uncertainty in slope (m/Am) is quoted to characterize the correlation between equivalent width and continuum 
flux. Also given are probabilities of the null hypothesis (P) from the Spearman rank correlation test. 
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Fig. 2.— (2c) 
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-5000 5000 

Velocity (km s" 1 ) 

Fig. 3. — Separation of variable and non-variable emission-line components. (3a) The high-state (solid black line) and low-state emission- 
line profiles (grey) are shown for C IV after subtraction of the continuum flux. The low-state profile has been scaled by a factor of f .4 to 
match the BLR flux in the high velocity wings. At lower velocities the profiles diverge, interpreted as the signature of a non-varying NLR 
component. The solution to the NLR profile is plotted with a dotted line, along with a Gaussian fit (smooth black line), as described in the 
text. The intrinsic absorption components for the blue (b) and red (r) C IV doublet members are labeled; Galacitic lines are denoted with 
a "G". (3b) Similar analysis is shown for N V, Lya, and Si IV. The residual NLR profile for each ion is overlaid with a scaled template of 
the Gaussian fit to the C IV NLR profile for comparison. A detector artifact (art.) and the O IV] emission-line multiplet are identified on 
the Si IV spectrum. 




Velocity (km s ') 
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-2000-1500-1000 -500 -Q000 - 1500 - 1000 -500 -2000-1500-1000 -500 -0000 - 1500 - 1000 -500 



Velocity (km s~') Velocity (km a"') 



Fig. 4. — Normalized absorption profiles in mean low (left panels) and high (right) flux states. Unocculted flux levels (1 — C) from 
a model with an unocculted NLR, and continuum and BLR covering factor profiles derived from the Lyman lines are shown (thick grey 
lines in 4a and 4c; red in 4b). Spectral regions contaminated with other absorption or detector artifacts are plotted with dotted lines. A 
model of the strong Galactic H2 absorption in the FUSE spectrum is shown (blue line). Radial velocities of the three strong kinematic 
components are marked with dashed vertical lines. For doublets having both lines free of contamination with other absorption, the doublet 
solution to the covering factor is shown with a horizontal bar. For enhanced S/N, the C II spectrum shown is the average of the ground 
state (C II A1335) and fine-structure (C II* A1336) lines. 
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Fig. 5. — C III metastable level populations. (Top) The computed populations for the metastable 3 Pj levels of the 2s2p term of C+ 2 
are plotted as a function of electron density for temperatures of T = 16000 (blue), 20000 (black), and 40000 K (red). (Bottom) Population 
ratios for the J = 2:0 (solid lines) and J =1:2 levels (dashed lines) are plotted as a function of density for the same temperatures. 
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Fig. 6. — Absorption spectrum for C III* A1175 multiplet. The normalized mean STIS spectrum shows the individual multiplct lines are 
unblended allowing measurements of column densities for individual s Pj levels in the 2s2p term. The spectrum is plotted in the velocity 
frame of the multiplct line with shortest wavelength. The radial velocity of each transition in the multiplet corresponding to the centroid of 
the Si IV feature in the averaged STIS spectrum is marked with a tickmark above the spectrum, with relative oscillator strength denoted 
by line length, and the lower level of each transition is identified by the J value. A fit to the absorption complex is plotted with a dashed 
line, with best-fit column densities for each level printed below the spectrum. 
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Fig. 7. — Photoionization modeling results for each intrinsic absorption component. (7a) For component la, solutions in U, Nh 
corresponding to the measured column densities in the low-state are shown as contours for each ion. The widths span the range of solutions 
based on estimated uncertainties. Regions of overlap give models that simultaneously fit those ions. The high-state solution to Si IV is 
included (unfilled lines) by shifting the solution along the x-axis by log(£/)=— 0.23, corresponding to the continuum flux ratio of the high 
and low states. This allows a direct comparison of solutions in different flux states. Model solutions for the measured limits (plotted with 
single lines) for component lb are shown on the right. Solution for components 2 and 3 are shown in Figure 7b and 7c. 
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Fig. 8. — Timc-dcpcndcnt population computed for Si IV in component 1 for various amplitudes of flux variations. Results for increases 
and decreases in ionizing flux by factors 1.4, 1.7, and 3 are shown. Flux changes were assumed instantaneous at t=0. Initial conditions for 
the calculations are from the best-fit Cloudy model for this absorber, with the density derived from the C III* A1175 feature. 
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Fig. 9. — Computed variations in the component la Si IV column density during the intensive monitoring observations. (Top) Results arc 
shown for two initial low-state models: the best-fit solution based on simultaneously matching the low and high-state measurements (solid 
line) and a higher-ionization solution from Figure 7a that matched only the low-state column densities (dashed line). The time-dependent 
populations were calculated between each epoch assuming step-function variations in ionizing flux equal to the observed UV continuum 
variability. Measured column densities are plotted with circles and error bars. The overall timescales and magnitudes of variations are 
matched well by the low-ionization model, but the variability is severely overestimated by the high-ionization model. (Bottom) The UV 
continuum light curve is shown for comparison. 
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Fig. 10. — Time-dcpcndcnt C IV and N V populations computed for component 2. Initial populations are from the best-fit low-state 
model (see text), with a step function increase xl.7 in ionizing flux at t =0. Populations for C IV (solid curve) and N V (dashed curve) 
are given relative to initial values, plotted as a function of the product of time (days) and electron density (10 4 cm -3 ). Lower horizontal 
lines denote the equilibrium values for the high-state models. Upper horizontal lines mark the 50% variability levels, which arc used to 
determine lower limits on the density based on the time separation between low and high states. 



